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Surface plasmon on graphene at finite T 
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Using the recently developed polarization tensor in (2+1) dimensions for the electronic excitations 
of graphene, we investigate the influence of temperature on the surface plasmons on graphene. We 
consider non-zero mass gap, but zero chemical potential. Plasmons may exist for both polarizations, 
TE and TM, of the electromagnetic field. For TE, the momentum region, where the dispersion 
function is real, appears bounded from below, whereas for TM it is bounded from above. We 
discuss the similarities of these features with those found previously in the case with non-zero 
chemical potential, but zero temperature. 


I. INTRODUCTION 


Electromagnetic waves travelling along a surface or an 
interface between different media, or along a thin sheet, 
are a well known phenomenon, especially on the surface 
of metals [1]. The amplitude of such waves decreases ex¬ 
ponentially in the directions perpendicular to the inter¬ 
face. Polarization and dispersion of such waves depend 
on the properties of the medium behind the interface or of 
the sheet. Such waves are frequently called surface plas¬ 
mons. However, there is no unique terminology. So, for 
instance, these waves are called sometimes simply ’plas¬ 
mons’. This latter notion is mostly used for waves of den¬ 
sity oscillations within a medium, a plasma half-room for 
example, together with their accompanying electromag¬ 
netic helds. These are longitudinal, whereas the surface 
plasmons are transversal. 

Surface plasmons may travel on graphene. These have 
been investigated since [2, 3] in quite a number of pa¬ 
pers [4-8] and in [9, 10]. The electronic excitations in 
graphene are modeled by a Dirac equation. The basic 
method of investigation is the random phase approxi¬ 
mation for the density-density correlation function. An¬ 
other, equivalent but possible less intuitive, method is the 
application of the methods known from relativistic quan¬ 
tum held theory. It rests on the description of the system, 
consisting of the electromagnetic held and the electronic 
excitations, as a quantum electrodynamics with photon 
in (3+1) dimensions and a fermion in (2+1) dimensions 
(with different speeds of light). Within this formalism, 
the problem to be solved consists of the calculation of the 
photon polarization tensor as a one-loop Feyn¬ 

man diagram using standard methods. The fermion has 
the Fermi speed v ^ c/300 as its ’speed of light’. This 
way the theory has two Lorentz groups, one in (3+1) di¬ 
mensions with the speed c for the photon, and another 
one in (2+1) dimensions with speed v for the spinor. 
These two interact through the usual minimal coupling 
in a consistent way. 

The polarization tensor n^^(/c) in the form, allow¬ 


ing an analytic continuation to the real frequency axis, 
was applied to the calculation of optical properties of 
graphene in [11], and to the calculation of the Casimir 
effect in [12-14]. For the case of graphene in a mag¬ 
netic field, see [15] and references therein. The relation 
of these calculations to those using density-density cor¬ 
relation functions, was discussed in [16]. A characteristic 
feature to be mentioned is that the polarization tensor 
has, in the presence of a medium, i.e., with T and/or /i, 
two independent form factors. The results obtained us¬ 
ing both methods and their relation to the van der Waals 
force were compared in [17]. 

The polarization tensor with both form factors was (re- 
)derived in [10], for ^ 0, T = 0, where also its relation 

to the reflection coefficients on graphene was discussed in 
detail. These formulas were applied to the investigation 
of surface plasmons. The case with /r = 0, T 0 was 
considered in [18] with application to the large thermal 
contributions to the Casimir force on graphene. 

In the next section we introduce the necessary nota¬ 
tions. In section 3 we investigate the TE-plasmon, in 
section 4 the TM-plasmons and conclusions are given in 
the last section. 

Starting from here, throughout the paper we use units 
with a = e^/(47r) for the coupling and v = vfIc for the 
Fermi speed, and refer to a = 1/137 and v = 1/300 as 
the physical parameters. 


II. POLARIZATION TENSOR IN PLASMONIC 
REGION 


Our notations follow [10], where electrodynamics with 
polarization tensor was considered in detail. Reflection 
and transmission coefficients for the scattering of a inci¬ 
dent plane electromagnetic wave off a graphene located 
at z = 0 are 
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correspondingly with TX—^-TE for the transverse electric 
(TE) polarization and with TX—J-TM for the transverse 
magnetic (TM) polarization of the wave. In terms of the 
form factors we get 

Qte = + .3. 

Qtm = 2^noo(i5), 

where ntr(p) = (p) is the trace of the polarization 

tensor. 

The photon wave function has plane waves propor¬ 
tional to ~ exp(—iwt + zkx ± ipz) (outside the plane, 
i.e., for z ^ 0) with k = {ki,k 2 ) (and k = |k|) - the 
momentum components parallel to the plane at z = 0 , 
and 

p = \/u;2 -k2, (3) 

for the momentum perpendicular to the plane. The ar¬ 
gument of the form factors is the momentum 

p = \JlzP- - {vkY, (4) 

which involves the Fermi speed v ~ 1/300. 

Surface plasmons appear as poles of the reflection co- 
efhcients, thus as solutions of the equations 

1 + Qtx = 0. (5) 

These solutions determine the frequency of the plasmons 
as function of the momentum, 

w = Wsf(k), ( 6 ) 

i.e., the dispersion of the surface plasmons. In order not 
to decay in time, these solutions must have real w. Fur¬ 
thermore, these must have imaginary p to decay in the 
direction perpendicular to the plane, 

p = irj = i\/ -I- (^ 7 ) 

This implies a dispersion below the continuous spectrum, 

vk < oj < k. ( 8 ) 


The lower bound comes from the polarization tensor, 
which otherwise would not allow for solutions of (5) with 
real w. In this way, it is just the second speed of light 
which opens a window for the existence of surface plas¬ 
mons. We call this region the plasmonic region. 

We represent the quantities Qtx, entering (1), in the 
form 

1 f p^ \ 

QtE = ^ I 2^2 (^00 90o) + {htr + fftr) 1 , 

QtM = 2y2J^2 (^00 30o), ( 9 ) 

where 

hoo = -a(f)%(p), (10) 

v = 2 a $(p). 
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FIG. 1. The plasmonic region defined hy vk < oj < 
Min(fc,cJs). The curve wj! is shown for g = 1, m = 1, 
V = 1/30. 


are the zero temperature parts. These depend on a single 
function, 

<I>(p) ( 11 ) 

2mp — {jp + 4 TO^)arctanh 2 ^^ , (w < Wg) 

2mp — {p^ + 4m^) ^arctanh^^ -b > ^s) 

only. This function has a threshold at 

uJs = a/ (vk)^ + ( 2 m) 2 . ( 12 ) 

Above, i.e., for oj > cOg, electron-hole creation is possible 
and the polarization tensor has an imaginary part. The 
resulting plasmons would not be stable. This way, 

OJ < ojg, (13) 



poses another condition on the plasmonic region. 

The temperature dependent part in (9) was derived in 
[ 11 ]. Changing notations to those used in this paper, it 
is given by an integral. 




dq 


1 


F e^/T’ -b 1 



(14) 


where ’A’ stands either for ’00’ or for ’tr’. Further nota¬ 
tions are 

F = Q _ 2vkq, 

Q =-p'^ - 2 XojT, N = sign((5)\/Q2 - a^. 

Moo = 4F2 +f+ 4AwF, Mtr =f+ 4m^ (15) 

These functions are real in the plasmonic region. 

It should be mentioned, that the denominator, N, is 
real and non-zero in the plasmonic region. In fact, it 
is real in a larger region defined by w < where 
= yjk'^v'^ — 2kqv + Tm? + q^ + \Jm? -b q^ {q is the in¬ 
tegration variable in (14)) is a solution of fV = 0, and 
uj < wj” holds. The plasmonic region is shown in Fig. 1. 
It does not depend on a. 

In this way, all quantities entering Eqs. (5), are de¬ 
fined, and in the next sections we investigate their solu¬ 
tions numerically. 
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FIG. 2. The starting point ko of the TE surface plasmon as 
a function of temperature T for a — 1/137, v = 1/300. The 
dashed line is at fc = 2m/-v/l — 


III. THE TE PLASMON 

We start with the investigation of the TE plasmon. At 
zero temperature the TE plasmon solution was discussed 
in [9]. It exists for all k. At finite T, this changes and the 
dispersion has a starting point, fep. It is shown in Eig. 2 
as a function of temperature. Here T and k are given 
in the units of mass. Eor example, at m = 0.025 eV we 
have T = 1, which corresponds to room temperature. 

The starting point of the dispersion curve is on the bor¬ 
der of the continuous spectrum, w = fc, and the endpoint 
approaches uj = ojs- Fig. 3(left) shows the dispersion rela¬ 
tions for TE surface plasmons for different temperatures 
and with a = 1/137, v = 1/300. The solutions are rep¬ 
resented by nearly straight lines, whose lengths decrease 
with increasing temperature. Represented in the (fc,a;)- 
plane, they all stick together. Therefore we represented 
them in the (A:, u;/r)-plane. 

As mentioned, these lines become shorter with increas¬ 
ing temperature. There is a largest temperature when 
the length shrinks to zero. This happen when the begin¬ 
ning point reaches the intersection of the border of the 
continuous spectrum, u) = k, with the threshold, to = lOs- 
The numerical solution of Eq. (1), with oj = k and si¬ 
multaneously w = LOs, gives Tmax ~ 3.76. However, al¬ 
ready for T > 2.25 the TE plasmon is invisible on the 
plot. As known from T = 0, the curves have a knee for 
k > (see [9]). For coupling stronger than its 

physical value, the knee at the end of the curves becomes 
visible, see Fig. 3 (right). 


IV. THE TM PLASMON 

The TM plasmon is the solution of the equation 
1 + Qtm = 0, (5). It is absent for T = 0. The dis¬ 
persion curve starts at A: = 0. The curve has an endpoint 
whose position depends on the temperature T. It can be 


FIG. 3. TE surface plasmon frequency as a function of k/T. 
Top: The dispersion of TE plasmon for physical values of 
the parameters and for several values of the temperature, 
T = 0.1, 0.3, 0.5, 0.7 (from bottom to top). Bottom: TE plas¬ 
mon frequency at T = 0.5 for different values of the coupling 
constant, a = 1/137, 1/30, 1/10, 1/5,1/2, 1 (from top to bot¬ 
tom). 


calculated analytically. 



(16) 


with 

_2 t ^2 o ^^2 

:F{T) = - -+-PolyLog[2,-e™/^]+Tln[l + e™/^]. 

6m m 

( 17 ) 

At zero temperature J^{T) vanishes, the endpoint is on 
k = 0 and the TM plasmon disappears. The dispersion 
of the TM plasmon is shown on Fig. 4. 

With increasing T the curves may partially go be¬ 
yond the region allowed for real solutions. Thus, two 
real branches of one solution appear, separated by a gap 
where the latter is not real. At relatively high temper¬ 
ature (T > 2.75), the real solution disappears at the 
boundary uj = ujs- It re-appears as second branch at 
large momenta, k < kend- For T = 2.75, kend/ni « 3176, 
which is far above the applicability region of the Dirac 
model. In this sense, for high temperature only the low 
momentum branch terminating at the threshold w = Wg, 
Fig. 4 (right), is physical. 

V. CONCLUSION 

In the foregoing sections we investigated the surface 
plasmons on graphene at finite temperature and non¬ 
zero mass (gap parameter). A common feature is that 
these plasmons may exist for both polarizations, TE and 


CJ 
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FIG. 4. The dispersion of the TM plasmon for physical values 
of the parameters and m = 1 for several values of the temper¬ 
ature, low T in the upper panel, higher T in the lower panel. 
The dashed lines are uj = vk and cu — k, the dotted line is the 
threshold uj = oJs- 


TM. Examples are shown in Figs.3 and 4. Thereby, with 
increasing temperature, the TE dispersion curve has a 
starting point ko, see Fig. 2, whereas the TM dispersion 
has an endpoint, see Fig. 5. As a result, for TE there 
exists a temperature, above which there is no solution, 
whereas in the TM case the solution disappears for low 
temperature. 

It is interesting to observe similarities with the case of 
chemical potential /x, which was considered in [10] at T = 
0. First of all, the plots for the dispersions look similar, 
see for TE Fig. 5 in [10] and Fig. 3 here, and for TM Fig. 
3 in [10] and Fig. 4 here. Also, for TE we have in both 
cases a starting point, see Fig. 4 in [10] and Fig. 2 here. 
In TM polarization there are endpoints in both cases, 
which have even analytic expressions, eq. (59) in [10] and 
(16) above. The corresponding plots are in Fig. 5. These 


show a similar behavior, i.e., starting in fcend = 0 and 
growing up with the corresponding parameter, /x resp. 
T.It would be interesting to investigate this similarity in 
more detail by considering surface plasmons with non¬ 
zero both, T and /x. 

Another interesting point would be the extension of 
the above methods to complex solutions, i.e., to reso¬ 
nances which can be viewed as plasmons with finite life 
time. These would be extensions of the dispersion curves 
beyond their endpoints. We mention that the extension 
beyond the starting points would give a solution in the 
continuous spectrum. These would have a wave function 
not decreasing in the direction perpendicular to the sheet 
and would be at best weakly localized at the sheet. 

^end 

m 



m 



FIG. 5. The endpoints of the dispersion corves of the TM 
plasmons for physical values of the parameters. In the upper 
panel for T = 0 /x 7 ^ 0 and in the lower panel for T 7 ^ 0 and 
/X = 0. 

Another application could be the question on whether 
the plasmons on graphene will play the same role as with 
plasma sheets with respect to the van der Waals forces 
between two sheets which was investigated in [19] and 
[ 20 ]. 
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